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ABSTRACT 


The Delaware Method of rating baffled shell and 
tube heat exchangers is used as the basis for a simplified 
mathematical model for predicting shell side heat transfer 
coefficients and pressure drops. Multiple non-linear 
regression analysis 1s used to express the heat transfer 
coefficient and pressure drop as a product of terms 
involving the five basic geometric parameters of a shell- 
and-tube heat exchanger, over a suitable range of these 
variables. 


A brief sensitivity analysis of heat transfer per 
unit pressure drop is presented, along with a proposed 
routine for determining exchanger size when the shell- 
Side coefficient is the controlling factor. A comparison 
is made between this method and the Delaware Method. 


Thesis Supervisor: Warren M. Rohsenow 
Title: Professor of Mechanical Engineering 
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NOMENCLATURE 


heat transfer area for tube side, (si) 


heat transfer area for shell side, (ft*) 


Speciite: nedtee: shelivmside fluid, (Btu/lbm - OF) 


(in) 


Shell inside diameter, 


diameter of outer limit of tube bundle, 


(in) 


tube outside diameter, 


(in) 


tube inside diameter, 


(in) 
EraeerOonwoOnreetoes that are in crosstflow 


fraction of total crossflow area, available for 


bypass flow around tube bundle 


Caatiedttomaletretion factor for flow across 
an ideal tube bank 


gravitational, conversion GOncmales ss. il) “x 10° 


lbm-ft/lbf-hr 


tube side heat transfer coefficient (Btu/hr-£t’-°F) 


ideal shell side heat transfer coefficient, 
(Btu/hr-£ft4-°F) 


exchanger shell side heat transfer coefficient, 
(Btu/hr-ft*-°F) 





AT) 


Shell side heat transfer coefficient correction 
factor for bundle bypass effects 


shell side heat transfer coefficient correction 
factor for baffle configuration effects 


shell side heat transfer coefficient correction 
factor for baffle leakage effects 


Shell side heat transfer coefficient correction 
factor for adverse temperature gradients 


icimciiemicaitaryanonen COCiLIlelent correction 
factor for the effect of unequal baffle spacing 


Colburn j-factor for ideal tube bank 


thermal conductivity of tube side fluid, 
(Btu/hr-ft-9°F) 


CieGHialeconatecrelvity Of Shell side fluid, ~ 
(Btu/hr-ft-9°F) 


thermal conductivity of tube wall, (Btu/hr-ft-°F) 
logarithmic mean temperature difference 
effective tube length between tube sheets, -(ft) 


baffle cut distance from baffle to shell inside 
diameter, (in) 


baietlter Spacine, (in) 





baffle spacing at inlet and outlet 
respectively 


of exchanger, 


number of baffles in exchanger 


number of tube rows crossed during 
one crosstlow section 


flow through 


number of effective crossflow rows in each 


window section 
total number of tubes in exchanger 


pressure drop for flow across one ideal crossflow 
section, (1bf/in24) 


pressure drop for flow through one ideal window 
section, (lbf/in2) 


Bue p heen, | 1m) 


tube pitch normal to flow, 


Gin) 


tube pitch parallel to flow, (in) 


total heat transferred in heat exchanger, 
(Bieuy ar) 


MRessureadrep COLEection factor for effect.of 


bundle bypass 


MucssMmcmanrOpMGcOrrectlon fLactor for effect of 


baffle leakage 


PORserrect Of 
Spacing 


HaAeceor 
baffle 


pressure drop correction 
unequal inlet and outlet 





Re shell side Reynolds number 
Ss 
S crossflow arga at centerline for one crossflow 
fs Seet lone. J( in) 
Sh shell-to-baffle leakage area for one baffle, (in*) 
Seb tube-to-baffle leakage area for one baffle, (in) 
S, area for flow through window, (in*) 
Swe total baffle window area, (in*) 
St window area occupied by tubes, (in*) 
U overall heat transfer coefficient base on shell 
° side heat transfer area, (Btu/hr-ft*-°F) 
We mass flow rate of shell side fluid, (l1bm/hr) 
W.. mass flow rate of tube side fluid, (lbm/hr) 
6 diametral clearance between shell and baffle, 
sb 
(in) 
van diametral clearance between tube and baffle, 
(in) 
Sctg approximate diametral clearance between shell 
anamolter tube dimit, (in) 
u shell side fluid viscosity at bulk temperature, 


(lbm/hr-ft) 
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Shell side fluid viscosity at surface temperature, 
(lbm/hr- ft) 


shell side fluid density, (lbm/ft”) 
baffle cut angle, (rad) 
Shell side Prandtl number 


tube side fluid viscosity at bulk temperature, 
(lbm/hr-ft) 
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I. INTRODUCTION 


A BACKGROUND 

The aim “of any heat exchanger rating method is to 
predict, as accurately as possible, the performance of a 
Meat exchanger with respect to total heat transferred and 
pressure drops encountered on both the shell and tube sides. 

The Delaware Method for Rating Baffled Shell and 
Tube Heat Exchangers (1) was developed from its inception 
in the late 1940's through the early 1960's as a result of 
numerous experiments conducted at the University of 
Delaware. This method is applicable to many variations 
on the baffled exchanger: multiple passes, finned tubes, 
various tube orientations, a wide range of Reynolds 
numbers, etc. As Bell has stated (1), the method does 
have its uncertainties, as all methods do. However, this 
method is widely used and easy to understand. 

Tinker (2) and Devore (3) propose similar methods 
fOr rating baffled exchangers, but, in the opinion of 
this author, are not as straightforward as the Delaware 
Method. Bell himself has proposed FOUcines fOr sizing 
baffled exchangers based on the Delaware Method 
(references 1 and 11), and this method has been used 
extensively in computer programs for "optimal'' exchanger 
determination. References (9) and (10) are two such 


examples. 
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The basic means of analysis in the Delaware Method 
takes into consideration five fluid streams and three 
heat transfer zones in the baffled exchanger, as illustrated 
in Figures 1 and 2. As can be’ seen by these simplified 
illustrations, the flow pattern for shell side fluid is 
complex at best. Reference (1) contains the complete set 
of equations for relating the heat transfer coefficient 
and pressure drops to the geometry of the exchanger. This 
reference also utilizes extensively graphs and tables 
for determination of some parameters. Mueller (7) has 
developed equations for many of these graphical relations. 
Appendix A is a compilation of the basic.equations of the 
Delaware Method and those derived by Mueller. This 


Sonprlation does have its limitations as noted in Section [.C. 





FIGURE 1 schematic Diagram of Streams on 
Shell Side of Baffled Heat 
ERenangers. 








PIGcURE Z(a) Region of Crossflow Between Baffle 
Virsmimternewcentral Bafile Spacing 





PUGURE Zi b) Flow Region Considered for Window Flow 





FIGURE 2(c) Flow Region for the End Baffle Spacings 
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Be PURPOSE AND METHODOLOGY 
PMicmuuUrpesesen tats Study 15 tO provide an insight 
into the relations between shell-side heat transfer 


coefficient and pressure drop and the five basic 


exchanger parameters: Tube diameter -- da; pitch-to- 
diameter ratio -- P/d.; bacrleseut ratio= - &/D;; batrie 
spacing ratio -- &_/D;; and shell diameter -- D,. 


It 1s anticipated that a simplified model for 
determining the heat transfer coefficient will be of 
considerable use to the designer as a "first cut" in 
Sizing an exchanger for a given job. Admittedly, this 
cannot replace the sophisticated computer programs 
available for the determination of "optimal" exchanger 
configuration for a given application. However, the 
method developed here will be a step in the right 
@enection for initial selection of parameters and will 
provide some "feel" for the effect of geometry variations 
on overall heat exchanger performance. 

This analysis was conducted utilizing the basic 
equations of the Delaware Method as modified by the 
assumptions and approximations enumerated in Section I.C. 
Utilizing these equations, those factors affecting heat 
transfer and pressure drop which are a function of 


exchanger geometry were calculated over a suitable range 





6 


of variables. The resulting values were correlated 
Mriiizing four FORTRAN subroutines (b): CORRE, ORDER, 
MINV and MULTR. The resulting non-linear correlation 
yields these factors expressed as a product of the 
five basic exchanger parameters, each with its own 
regression coefficient as an exponent, and a single 
correlation coefficient. These regressions provided 
acceptable correlation results, the maximum overall 


error being on the order of five percent. 


G. LIMITATIONS/ASSUMPT IONS/APPROXIMATIONS 

Ihe Limitations 

In this analysis, no attempt has been made to 
correlate data for the shell side heat transfer co- 
efficient where sealing strips are used to limit the 
bundle bypass flow. Nor has an attempt been made with 
regard to enhanced surfaces for heat transfer, i.e. 
finned tubes. Finally, only the single pass condition 
Bwmcic tube Side fluid 1S taken into consideration. 
Admittedly, these restrictions limit the scope of 
Enis study with regard to overall applicability to 
Meat exchanger design. However, it is hoped that this 
analysis will provide insight into the basic relations 
Between shell side heat transfer coefficient and the 


fundamental exchanger goemetric parameters. 
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oF Assumptions 


It is realized that manufacturing considerations 
Seeect Silgniticantily the desion, construction and final 
form of any heat exchanger. Reference (5) sets standards 
for the tube-to-baffle and baffle-to-shell diametral 
clearances for various applications. These values are 
predominantly of the "step function" type, in that they 
can not be expressed as a continuous dependent variable 
based on any of the independent exchanger parameters. 
For the purposes of this study, both clearances were 
assumed constant throughout the range of parameters 


considered. The values used for these clearances are: 


Shell-to-baffle clearance: 6 ob = 0.175 inches 


Tube-to-baffle clearance: Sy 0.03125 inches 


These assumptions again allow concentration on the effect 
of major exchanger parameters. 

The second major assumption is the neglect of 
muequal batfle spacing for inlet and outlet passages. 
This assumption is made in both heat transfer coefficient 
and pressure drop calculations. Unequal baffle spacing, 


in an exchanger where the ratios 
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contributes to an approximate 2-3% reduction in the over- 
all heat transfer coefficient for a sufficiently large 
number of baffles. Where pressure drop is concerned, 
the assumption must be made with the realization that 
the end zones be considered as crossflow areas identical 
to the zones between the baffles. Equations (19) and 
m5) of Appendix A reflect these changes. 

oF Approximations 

In order to implement any numerical analysis of 
the type used in this study, it 1S necessary to derive 
relations from data provided by the manufacturing 
community and from curves of friction factor and 
molburn j-factor. 

For the former, references (1), (5) and (8) were 
used for two such derivations. The number of tubes in a 
heat exchanger of particular shell diameter was corre- 
Maced aS a function of shell inside diameter, pitch- 
Meetube diameter, and tube diameter, for the single- 
pass condition. These results are: 


For square/rotated square tube layouts: 


1 le d. Z 
N. = 0.5559 ——_—_ _|— (1) 


2 
a 3 





LG = 


For triangular tube layouts: 


pn, 2-1797 74 )2 
(2) 


Mw = O.40"5e == — 

: a? P 
O 

These approximations give results within 5% over the range 
of shell diameters from 10 to 60 inches. The clearance. 
for bundle bypass flow (D. - Dog? was correlated as a 

function of shell diameter also, resulting in the 

approximation: 


0.3066 2 


§ = (D = (0.807 De (3) 


stz i ~ Poesy? 
Pictimesmonmon andes ane praphs of friction factor 
and Colburn j-factor versus shell-side Reynolds number 
for the various tube orientations at varying pitch-to- 
tube diameter ratios. Correlation results are presented 
maumetiaples 1 and 2 for friction factor and Colburn j-factor, 
respectively, for the separate ranges of Reynolds numbers 
shown. 
Table 3 presents relations between tube pitch, 
mabe pitch normal to fluid flow and tube pitch parallel 
to fluid flow for the various orientations. These 
relations are used in equations (1), (3), and (4) of 
Pependix A in calculating the various factors associated 


with heat transfer and pressure drop calculations. 
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P = tube pitch 
oe = tube pitch normal to flow 
Py = tube pitch parallel to flow 
Tube 
Orientation Pn Po 
square Ie IO) aie nD oe = 
rotated 
square ae 10] eee 
triangular P/2 73/2 P 


ABE. 5 Relations Between Tube Pitch 
P, and Pi» P tor Various Tube 


Ie 


Orientations 
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De NEGREGHREDTEACTORS Jy and J. 


For preliminary estimation and sizing as proposed 
by this simplification, the factors J, (correction due to 
adverse temperature gradients) and J. KGomrcctTonmeaic £6 
unequal inlet and outlet baffle spacing) have been 
mee lected. 


fieomacecoumtators the effect or J the final heat 


Ss 
exchanger length must be known, which allows determination 


of the end baffle spacing. 2 may then be calculated 


TOM ; 
2p Wes 49 aloes, 
Game) + 22s. + Q 
jo = wb seb 0 (4) 
(Ny -1) * aoa ay ~s 0 
where 
“ aa 7 2s alts» 50 ~ 5 al US 


ire initial value of ie must then be multiplied by Je to 
mye a more accurate estimate. 
iie ace on J. 1s calculated only in laminar flow 


(Re . 00) from: 





0.18 (5) 
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Should this condition (Re, CeiwOjmexiSt. again adjust A, 

By Multiplying by J. Ereounes Orand 7 show graphically 
* 

the relations for J, VS hei and J. VS No: These graphs 


may be used in lieu of actual calculations. 
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Geommcertenm aebor fOr eAaverse lemperature 
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BiweeeccnceollONSE ACTORS INFLUENCING THE 
Hene eeaNoceR CORREICIENT 


A. BUNDLE BYPASS FACTOR -- Jp 


The bundle bypass factor, Jy» accounts for the 
effect of fluid flow between the shell and outer tubes 
of the tube bundle. This portion of the flow bypasses 
the main crossflow section of the tube bank and has a 
Significant effect on heat transfer. Equations A-4, 
A-5, and A-18 are used to calculate Jy: where the shell- 
to-tube-bundle clearance (D; - Do+1) has been approxi- 
mated by equation (3). Equation A-18 was developed by 
Mueller (7) where for laminar/transition flow Con =) eo 
emcd for turbulent flow nn [ice wor DUrpOSeS Of this 
analysis only two regimes of Reynolds numbers are 
Semsitdered: 1 < Re < 10° fone laninanry/ transition f£1ow. 
and 10° < Re_ < 105 for turbulent flow. 

The non-linear regression results for Jy are 
presented in Table 4 and a representative curve is shown 
mirrFigure 8. As can be seen, Jy can be represented as 
a function of shell diameter D;, tube diameter da» and 
pitch-to-diameter ratio P/d.. 

From Table 4, it can be readily seen that tube 


diameter has little effect on Jy: Over the range of 


d. considered (0.5<d.<1.5), the change in Jy due to 
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increasing (or decreasing) d, PCO, tous 5 (or to 0.5) 
is less than 1%. 

Figure 8 is applicable only for turbulent flow 
Over square tube orientations. However, we may estimate 
Jy for other flow/orientation conditions easily from 
Table 4 and can estimate changes in Jy in the following 
fashion. If J, is desired for rotated-square layouts 


b 
in turbulent flow, we may form a ratio: 


Jy fc (p/a.) 2 (.) 2 
~  (C, (P/4,) ? ,) 32r 
a ; re. (Gavel "2 op "33, 
re tr/aj * tp,) 73s 
(6) 
Co s. - a oh, - *s. 
a an (P/d.) (D;) 
Og 


ims simple relation yields the result that rotated-square 
layouts, for P/d. = 1.5, D; = 30, result in a bundle 
bypass factor approximately 6% higher than square tube 
Orientations. 

Such insight into the intricacies of the Delaware 
fmaeneod Can allow the designer to make reasonable "fine- 


tunins” corrections to a preliminary design without the 


need for costly (and sometimes cumbersome) computer models. 
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B. BAFFLE CONFIGURATION FACTOR -- J. 
Foucemonom-ceamdeA-1G represent the batfle con- 
miguration factor, Jos which amounts for changes in heat 
transfer due to those tubes physically located in the 
baffle window and not subjected to full corssflow from 
the main stream of fluid. Equation A-16 is a linear 
Boproximation to the graph of Jo developed by Bell. Using 
Sauation (3) for Dot e Jo can be reasonably approximated 
as a product of shell diameter and baffle cut (or baffle 
cut ratio). This result is independent of flow condition 
and tube orientation. Regression analysis results in 
the following equation for Je! 


-0.5467 -0.0184 
On 2G? 2} (D;) (7) 


es 
Figure 9 is a family of curves plotting &/D; Vs J. at 
various values of D.. 

As is easily seen from this figure, the dependence 
ony J. on shell diameter is minimal when compared with 
Maat Of the baffle cut ratio. For the range of baffle 
cut ratios considered in this analysis (0.2 < &/D; <A 
the heat transfer coefficient may vary from a 15% increase 
to a 25% decrease. Bell (1) states that for "well-designed" 
exchangers, a Should be about 1.0, which corresponds to 


a value of &/D; Omer pOwe se e)2 5 . 
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It 1s interesting to note that in the optimization 
routines of references (9) and (10), the baffle cut ratios 
for "optimal"' design were 0.25 and 0.2, respectively. 
Admittedly, the objective functions for these optimizations 
were Slightly different (heat transfer area and total 


tube length, resp.), but the similarity of results is 


E1onificant. 
c. Poapeie benkAGE FACTOR -- J% 


The baffle leakage fa Geom . is by far the most 


Q? 
complex of all the factors affecting the overall heat 
meransrer coetficient. It accounts for the effect of 
leakage through the baffle due to clearances between 
baffle and shell and between tube and baffle. In this 
memard, lt 18 a function of all five basic exchanger 
parameters. Equations A-6, 7, 8, 15, 16, 17 must all 
be used in the calculation of Jo; where some additional 
dependence is noted as a function of tube orientation. 
The non-linear regression results for Jy are 
presented in Table S. Numerous graphs and families of 
curves can be generated using these results. However, 
for the sake of simplicity, Figures 10a, b, and c are 
represenative of the effect on Jy of several parameter 


Variations. Bell states that the value of Jy should 
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FIGURE 10b Baffle Leakage Factor J, vs P/d 
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FIGURE 10c Baffle Leakage Factor Jy VS P/d. 
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be 0.7 -- 0.8 for weli-designed exchangers. From 
Figure 10, one can easily read an appropriate value 
of the necessary parameter that will be a reasonable 
meartiioe POlnt Lor heat exchanger design. For other 
Shell diameters and tube diameters, the relations of 
Table 5 may be suitably manipulated to arrive at 


reasonable values of the necessary parameters. 


Di MUO mE mCOnkE EAL LON RESULTS 

To facilitate the overall prediction of the 
Meat transfer coefficient as a function of those 
variables considered, the actual shell side coefficient, 
h 


was calculated as a product of h. and the 


om ideal 


aforementioned factors: Jy Jo» Jy. Colburn | —Lactor 
correlations of Table 1 were used for the various flow 
regimes as were correlations for N. and 5 ste: The 
fluid properties and mass flow rate terms were factored 
out of the equation and only those terms involving 
geometric parameters were used for the overall cor- 
relation. These terms were then correlated using the 
aforementioned Fortran subroutines. The correlation 
results are presented in Table 6 for the various flow 


regimes and tube orientations. It should be noted that 


all input dimensions for the exchanger parameters are 





=47= 


in inches and that the correction factors inherent in 
these equations have been accounted for in the cor- 
memration constant. The resulting units for Ay are 
those typically expected of the English system, Btu/HR- 


rT? .° 


F. The units for fluid properties are as listed 
im the nomenclature. 

Table 6 is a "busy" listing and many graphical 
representations can be drawn plotting the heat transfer 
coefficient A, against the various parameters. That 
Will not be done here but some aspects of Table 6 are 
noteworthy. 

Most exchangers are designed for operation in 
turbulent flow and for this regime it can be seen that 
; A, is nearly inversely proportional to shell diameter, 
the exponent ac being close to -1.0. This is a striking 
result, in that for desired increases in Ao» Changes in 
shell diameter maybe the most rapid means of attaining 
tne desired value. 

The dependence of A, Onebatile Cut ratio, 

2 /D;, is independent of flow but dependent on tube 
Orlentation, and Vs Varies nearly as the square-root 
of & /D;. Pitch-to-diameter ratio, P/d., has little 
effect on square tube orientations, but can cause 
Significant changes in triangular and rotated square 


patterns. The dependence of Ae On tube diameter and 
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battle Spacing ratio 1S Straightforward and needs no 
further explanation. 

PROMelamke G6 ot 1S quite easy to predict the 
heat transfer coefficient given the exchanger geometric 
parameters and fluid properties. Only a quick Reynolds 
number calculation is necessary beforehand to determine 
the applicable flow regime. Changes in Ay are pre- 
dicted simply by calculating the ratio of the parameter 
varied raised to the appropriate exponent. This is 
quickly accomplished if only one parameter is varied 
Mera @1iven time. For large changes one must be cautious 
that the flow regime remains as desired, i.e., turbulent 


mae laminar/transition. 
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PUeomece Ono aerecl ING PRESSURE DROP 


A. BAFFLE LEAKAGE FACTOR -- Rx 

The pressure drop correction factor for baffle 
leakage, R,, accounts for pressure drop changes in much 
the same fashion that Jy accounts for changes in heat 
transfer. It relates the leakage through the baffle due 
to tube-to-baffle and baffle-to-shell clearances to the 
geometric parameters of the heat exchanger. Equations 
fee oe OO. 7, S, 145 15, 20 and 21 are all used to 
eaEculate Ro: Dependence on tube orientation is accounted 
for in equations A-4 and A-6 by the use of the correla- 
—E1ons for N.. iieequartons (1 )\) ana (2) and the tube pitch 
imenhations of Table 3. 

The non-linear correlation results are presented 
in Table 7 and a representative graph as Figure ll. 
Baeures 10 and 11 show a striking resemblance of forn, 
which is to be expected. One can logically assume that 
an increase in heat transfer coefficient due to less 
leakage through the baffle will also result in increased 


Deessure drops. 
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les BUNDLE BYPASS FACTOR -- Ry 

The bundle bypass factor, Ry» accounts for pressure 
drop corrections in the same manner that Jy accounts for 
heat transfer changes. From Equation A-22, it can be 
seen to relate that correction to Fsbp’ the fractional 
crossflow area available for bundle bypass. 

Table 8 presents the correlation results for Ry: 
Similar to Jy is the weak dependence on tube diameter 
and the strong dependence on pitch-to-diameter ratio. 
This comes as no surprise since increases in P/d. "open 
up'’ area for flow through the crossflow bank, thereby 
decreasing the fraction of bypass area. The increase 
in Ry with increasing D; is less than that of P/d., 
Since the shell-to-outer-tube-limit clearance (Equation 3) 
also increases with D.. 

Figure 12 is a graphical representation of Ry 
vs . ee at various D:, in the same manner as Figure 8. 
Meesnould be noted that over the range of qd. considered 
io. < d, << leopes thie value of Ry will decrease about 
5% with increasing d.> a much larger percentage decrease 
than can be predicted for Jy (and subsequently, heat 


transfer). 
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c. CROSSFLOW AREA PRESSURE DROP 

Equation A-25 gives the ideal pressure drop for 
flow across the tube bank between baffles. When 
multiplied by Ry and Ry» the effects of bundle bypass 
and baffle leakage are accounted for. This quantity 
has been correlated in the same manner as the previous 
values. The friction factor correlations of Table 2 
were used along with other Ore ions for tube 
mecntation. Table 9 1s the resuit of these overall 
correlations for the two flow regimes and various 
orientations. Note that the exponent X¢ of the basic 
equation depends on flow and orientation as a result 
@reethe friction factor correlation. Hence, the 
Mmeessure drop can not be expressed only as a function 
of exchanger parameters. From Table 9 it is obvious 
mee alt the geometric parameters inversely affect 
pressure drop, a result that may have been intuitively 
obvious, but the extent of which can now be reasonably 
quantified. Essentially, the increase in any of the 
five basic parameters allows for more area for fluid 
flow with less restriction; hence, lower pressure 


drops. 
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DB . WINDOW AREA PRESSURE DROP 

Pudi etonme= 25 eaves the ideal pressure drop for 
the baffle window section of the heat exchanger. When 
moltiplied by Ro» the effect of baffle leakage is taken 
into account. Correlation results are presented in 
Table 10 for differeing tube orientations. Unlike 
crossflow pressure drop, this quantity is independent 
of flow regime. 

These results are very similar to the crossflow 
pressure drops in that pressure drop is inversely 
meopertional to all variables. However, in the window 
euea the relative importance of the parameters is slightly 
different. For the window area, shell diameter is the 
predominant influencing factor instead of pitch-to-diameter 
ratio. The relative influence of baffle cut and a baffle 
spacing ratios is nearly the same for both area pressure 
drops. Tube diameter impact is relatively the smallest 


for the turbulent crossflow area pressure drop. 
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IV. PERFORMANCE AND SIZING OF HEAT EXCHANGERS 

A. Soult) weANSL YSIS OF HEAT TRANSFER PER UNIT 

PRESSURE DROP 

For any given application, when determining the 
"optimal" size of a heat exchanger, some feel for the 
overall effect of parameter variations on heat exchanger 
Berrormance 1S needed in order to judge the direction 
in which the design variables should be altered to arrive 
at a reasonable result. The performance of a heat 
exchanger, or the criterion against which an exchanger. 
1s to be considered optimal, can vary widely depending 
on the depth of the analysis that is being performed. 
In most practical applications the ultimate criterion 
Mmeeeost. In this analysis, cost will not be considered 
but a factor which can influence cost to some extent 
will be analyzed. 

From the corrleation results for heat transfer 


coefficient, Table 6, and the basic heat balance: 


QO ae (8) 


consider the case where Uy 1S primarily dependent on Ao. 


Equation (8) then can be written as: 





Be 


Q 


At om HS, 


where Ky is defined in Table 6. But A = md Nip /12 (£07). 
and the results from Equation (1) (or (2)) can be used 


for N, to express A in terms of d., D,, P/d.. Deane, 


ya 
Q 
Oe as ————__—__—— 
AT om eal 
which yields: 
m™ do N. Ae 


This quantity is a non-dimensional expression for heat 
transfer per unit length of the heat exchanger. From 
the correlation results for crossflow and window area 


pressure drops, (Tables 9 and 10), we may forn, 


Similarly: 
“s 
r/L?-=-lh Oe Hee --.) 71. 10 
if a. (AP, 5 + OP, 5)/ (10) 


This 1s a simplified expression for the total shell-side 


pressure drop per unit length of heat exchanger. Then 





Rime 


t/t can be expressed as a quotient of those correlations 
already derived. It can be seen from all the correlation 
results that Me is not only dependent on the geometry 
of the exchanger but flow rate also, specifically, 

W/w,- For this analysis, we will consider only turbulent 
flow across square-oriented tube bundles. The applicable 
correlation results were selected, the ratio Th/Ty 

formed and that ratio was graphed as a function of the 
exchanger parameters. It should be noted that a Reynolds 
number of 10° Waseuisead. Gesulting im a value of 


5 


ey = 2.026 x 10° at the mid-range of the exchanger 


geometry variables (4, SO) a Bal Sed 35; a _/D; =i. Os 
Se= 55). EFieures 13, 14, 15, 16 and 17 are various graphs 


1 

ohe pl Oe From Figures iL gaily it 1s obvious that shell 
diameter D. is the most influential parameter affecting 
heat exchanger performance. Figure 13 graphically 
illustrates this influence for several values of baffle 
spacing ratio, &_/D;. This is not to say that we should 
always optimize an exchanger at large values of D.. 
Although r,/T increases, this graph does not allow 
meiieany CONStraints to be considered, ¢€.¢g. maximum she11- 
Side pressure drop. However, we may conclude that where 
large amounts of heat are to be transferred with low 


Shell-side pressure drops as a requirement, "large’ shell 


diameters will be required. 
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Figure 14 illustrates the influence of baffle 
Spacing on Tie eee Vemalyvalucs sor D.. This parameter 
also favorably affects the heat exchanger performance 
Bie certainly not to the extent of shell diameter. Again, 
the realization of the absence of a pressure drop con- 
straint must be accounted for. 

etic minitSenates. Gnesnelationship between 
Syl es and tube diameter. This is the only parameter (of 
those considered) which adversely affects Bad Baae The 
most basic explanation for this effect is that increases 
in tube diameter will more quickly decrease heat transfer 
than they will pressure drop. 

Figure 16 and 17 relate the effects of baffle 
eae ratio, &/D; and pitch-to-diameter ratio, P/d, to 
yl Bae It is obvious that at the conditions analyzed 
(turbulent flow, square tube orientation), neither para- 
meter has any significant effect on heat exchanger 
performance. 

The five figures shown can be of assistance 
when determining exchanger size for a given "job". The 
qualitative effect (or lack thereof) of the five para- 
meters will help guard against "false starts" in any 


iterative sizing scheme. 
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Iie ER@ueeeumol ZING KOUTINE FOR CASE OF SHELL-SIDE 
PRO@ieweres ONLY 


For the case where the shell-side heat transfer 
coefficient 1s one controlling factor (the tube side 
coefficient has little impact on overall Uo); it is 
@esirabie to determine the exchanger geometry for a 
feven “job': total Q transferred; AP. Speciticd, We» 
AT. (shell-side), and fluid properties all given. We 
may desire to "optimize" the heat exchanger on minimum 


volume where volume is simply expressed as: 





Oo 
ii 
pe 
& 
O 
4 


U. = h A= — a 


pacids 
ae ae (a2) 


OO 
il 
> 
J 


With knowledge of the fluid properties we may form the 


quantity rh: 


a (13) 
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where Ky Pmoeeeit loadin table oO, Correlation results 


for shell-side coefficient. Then 
T/L = 2 eee (14) 


We must now begin to specify the geometry of the 


exchanger. First, pick do: P/d., Ds as the beginning 


points. From equation A-12 for Reynolds number: 


ied. WwW 
R = Oo. 6S 
e 
S Ue Sn 
(15) 
W a | a a 
a. Ss 1 2 3 4 
eC ees / denen (25/0) ~ (D;) 
Ms 
where 
oa OZ. 224. a, = U7. a, = - oS. 
an rol eee, = er 9 35 
S 4 


we may either determine &/D; (if we desire to operate 
at a certain ie Ouaprek &_/D; and determine the flow 
regime (laminar/transition or turbulent). 

Once this condition has been finalized, the 


mreper correlation for ne is selected from Table 6, 
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Pmowine the selected tube orientation: The last para- 
meter to be determined is &/D;. As Bell has stated, 
for ''well-designed" exchangers, J. Sileiia oS elotonshe Ile 
which, from Figure 6, yields a value of &/D; of about 
0.25. The initial value of h, may now be calculated 
from the correlation results of Table 6 and exchanger 
length L determined from Equation 10, using the proper 
N. correlation for tube orientation. Volume may then 
be easily calculated. 

Pressure drop calculations can be made from the 
correlation results of Tables 10 and 11, knowing the 
exchanger length L. This must then be checked against 
the allowed AP .. iis AP . aS AP . , we may 

calc allowed 
decrease the size of the exchanger to utilize more of 
the allowed AP... This is most easily accomplished by 
decreasing D. and recalculating hy» lov AP . and V, to 
arrive at some minimum volume. 


If AP the procedure 1s 


Zicalie 7 OMS idenieal’ 
reversed, 1.e., increase D. and recalculate the 
necessary quantities until a convergence within some 
désired bound is obtained. As in the case in any 
Sizing procedure, some notes of caution and guidelines 


are in order. When increasing shell diameter (or any 


of the other parameters), caution must be taken that the 





oy, = 


flow regime remains as desired. It 1s also possible to 
Dhine tune’ the exchanger by altering the other geometric 
walues. If the difference NG <a OPS aeaed is small, 
we may desire to change the baffle cut or baffle spacing 
ratios. Small changes in these parameters will also 
result in changes to the shell side pressure drop and 
heat transfer coefficient, the effect of which can be 
examined by the exponents of the applicable correlation 
result. We may also use Bell's insight into the heat 
transfer corrections to assure that the heat exchanger 

1s reasonably well designed. For the heat transfer - 


eOorrection factors Jy» Ja» J Bedlecuyes wae following 


L, 


ranges as good guidelines: 


From Figure 9, 10, and 11, we can pick reasonable values 
ee &/D; and &_/D; as starting points for the exchanger 
geometry and, should changes be made, ensure that we 
remain in a normally expected range for those values. 

The use of those figures as aids to sizing should prevent 


haphazard variations and point to a reasonable solution 
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MiGUREw ES ortmpltired Logic Sequence for 
Proposed Sizing Routine 





SEG 


in an expedient manner. Figure 13 gives a simplified 
flow chart of the logical sequence for this proposed 


pizins routine. 


c. ExaMe les Or PROPOSED SIZING ROUTINE WITH 
' COMPARISON TO THE DELAWARE METHOD 


The preceding sections of this analysis have dealt 
With the actual results of the non-linear regression used 
to simplify the heat transfer coefficient and pressure 
m@eop calculations. In this section, a hypothetical 
example for determining exchanger size is given to show 
the actual mechanics of the proposed sizing routine and 
momcOMlpare those results with the calculations of the 
wenaware Method. This comparison is done for the case 
where only the shell-side heat transfer coefficient and 
pressure drop are of concern to the designer. 

The “job" considered in this example is as 
follows: 

il) Shell-side mass flow rate of water: 

iis = 250,000 lbm/hr 
2) heat water from 50°F to 90°F 
5)) Maximum shell-side pressure drop: 


AP. < 8.0 psi 
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The following fluid properties are assumed: 
1) ¢. = 1.0 btu/1bm - OF 
2) cee Peo oom, att. | ac sale ds sumed 


4 = 100°F) 


wall 
: ; sere 
Sa) Me eee oem, Mitt | at aan LOE} 


S 

4) Bape ee US 

5) p. = 62.075 1bm/ft” 

We desire to operate in the turbulent flow regime 
and must check the Reynolds number for our exchanger. 
ime kilg tube diameter d. ="O0R75 in, pitch-to-diameter 
ratio P/d. less battle Spacing ratio &/D; = 0.75 
and shell diameter D. =o wml a wcneck Of Reynolds 
number using Equation (15) yields Re. Se oo soow Nebel 
ensures turbulence. Using a square tube orientation, 


baffle cut ratio & /D; = 0.25, we select the applicable 


@errelation result from Table 6: 


h, We -. 355 Me Oils 
—=w, |= (4.07) (d,) (P/d.) 
Ky Us 
-, 4492 - 376 -1.032 
(2./D;) (2./D;) (D;) 
where ; 
= ae hea 
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Using the values given results in a heat transfer 


coefficient of: 


0 


hee otos. 5 Ben/hr-£t Ig 


O 


From Q = AU, AT. - Ah AT, o> where 


eee ands N. = 0.3559 = 
OQ 


ra, N, L 5 EoLTSAE pa He 
12 


Mieelengtn Lb of the heat exchanger is calculated: 


ee = 9620.99) “EF e:, 
then volume 


Var eS ete | 
Using the correlation results for crossflow and window 
area pressure drops from Tables 9 and 10 respectively, 
yields the total shell-side pressure drop: 


AP . = .29 psi 


AS a comparison, utilizing the Delaware Method (Reference 


1), the following is presented: 





Oe 


Proposed Method 


Lae 


ha CStu/hr-tt)- F ) 463.5 
mec rt ) 6.99 
5 = 
eeceit™ ) 218 3 
AP_. (psi) OR aS, 


Delaware Method 


454.8 


WAS 


Zan OS 


Although the results are quite close, we are far from 


"optimized" in using the allowed shell-side pressure drop. 


We may decrease shell diameter D. to most easily 


eecomplished this. For the next interation, pick 


D. = 15.25 in. This "odd'' size is picked only because 


Bell's tables in Reference (1) list all the applicable 


parameters for this size (Doty t? 


N etc.) and inter- 


polation between those listed values would not result in 


a meaningful comparison. 


A recalculation of the heat exchanger using this 


new value of D. yields the following result: 


Proposed Method 


ho (35% 16 
L 12.09 
V LS eo 


AP oO 
S 


Delaware Method 


789.0 


Ie Oe 
ners) S 


Sc (eat 
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These results are as agreeable as those of the 
first iteration, yet the exchanger may be further de- 


Breasea in size. A third iteration, using D. Sr Sac or 11) 


fave Las: 
ERowosed ietnioed Delaware Method 
A, 907.8 O17.5 
i, 14.20 3 14.61 
V IES 13.99 
AP . Tx 55 irs | 


These results show that we have come close to a 
preliminary size for the heat exchanger that meets the 
required "job". Some comments concerning the previous 
calculations are noteworthy. 

In decreasing shell diameter in these iterations, 
Mmece recalculation of A, by the proposed method is 
accomplished simply by multiplying the previous value 
by the ratio of shell diameters raised to the applicable 
@Momrelation exponent. Pressure drop calculations are 
estimated in much the same manner, except that the 
crossflow and window area terms must be done separately. 
ime only ‘lengthy’ calculation involves that of exchanger 


Hength, L. 
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For the sizing determination by the Delaware 
Method, no less than twenty-one (21) separate calculations 
and seven (7) graphical readings must be done for each 
meeration. Gor this proposed method only five (5) 
calculations were necessary: heat transfer coefficient 
(Table 6), overall length (Equation 12), pressure drops 
for crossflow and window sections (Tables 9 and 10), 
and total pressure drop (Equation A-25). 

Mimemotrek, handy method 1s, in the opinion of 
tars author, much more preferable to the original method 
for preliminary heat exchanger sizing and "back-of-the- 
envelope" calculations. It is realized that no simplified 
heat exahanger model can replace the more sophisticated 
computer optimization methods that have been and are 
being developed. However, in a relatively short time, 
the designer can achieve a realistic estimate of the 
heat exchanger performance and size, and can use these 
results as possible inputs to the more elaborate routines. 
The closer an initial estimate is to the actual optimun, 
the more time and money can be saved in the more 
sophisticated (and expensive) methods. Minimization 


of cost will almost always be a major goal. 
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APPENDIX A 


DacSmemeciiatrols USed Lon heat transfer coefficient 


ama pressure drop calculations. 
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2 Ge 
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(16) See 72 F 
7 ) Jo = 0.44 (1 - ee) ee eee clea ee rjle 
(18) J, oo “bh Fsbp- 

(19) h =h Sees aie 


(20) DS = (5st (hes oy ee Wat 


Do ees Ee p 
(21) a a LoSS CL Sard yict a 


Ps 
(22) 2 2 & 392 





=7-6=- 


(24) AP . = 


(25) AP, = CCN, + 1) AP, G RL +N, APY GIR, 
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